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Abstract. We present an analysis of the space motions of 742 subdwarf stars based on the sample of Carney et al. i!994l . 
Hipparcos parallaxes, TYC2+HIP proper motions and Tycho2 proper motions were combined with radial velocities and metal- 
licities from CLLA. The kinematical behavior is discussed in particular in relation to their metallicities. The majority of these 
sample stars have metal abundances of [Fe/H] > -1 and represent the thick disk population. The halo component, with [Fe/H] 
< -1.6, is characterized by a low mean rotation velocity and a radially elongated velocity ellipsoid. In the intermediate metal- 
licity range (—1.6 < [Fe/H] < -1), we find a significant number of subdwarfs with disklike kinematics. We interpret this 
population of stars as a metal-weak thick disk (MWTD) population. 
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1. Introduction 

Studies of the kinematics of various stellar populations in the 
Galaxy, in particular the thick disk and the halo, have long been 
limited by the availability of large samples of stars with mea- 
surements of proper motions, radial velocities, distances, and 
metallicities. Such data are required in order to constrain plau- 
sible scenarios for the formation and evolution of the Milky 
Way. Samples of nearby subdwarf stars with high space mo- 
tions provide an observationally convenient probe of the struc- 
ture of the Galaxy. The large, proper-motion-selected, stellar 
samples of Carney, Latham, Laird, and Aguilar ( 1994 here- 
after CLLA) have proved particularly valuable for studying the 
kinematics and chemical abundances within a few kiloparsecs 
of the Sun. 

The correlation between kinematics and metallicity gives 
useful information for formulating theories of galactic struc- 
ture. Differences in chemistry and space velocities are crucial 
in defining the different populations within the Galaxy and in- 
ferring their origins. Relevant studies of the kinematical behav- 
ior of stars ,in particular in relation to their metallicities ,were 
presented by e.g. Morrison, Flynn, & Freeman ( 1990 hereafter 
MFF) using a sample of K giants whose metallicities are mea- 
sured using the DDO photometric system, Nissen & Schuster 
dl991> using late F and G dwarfs and subgiants, Chiba & Yoshii 
)8 1 using red giants and RR Lyrae stars, Martin & Morrison 



1998 1 with a sample of nearby RR Lyrae stars and Chiba & 



Beers (2000) using 1203 metal-poor solar-neighborhood stars. 

The Galactic halo is characterized by a roughly spherical 
space distribution with close to zero net rotation. Its stars are 
metal poor, with a peak metallicity at [Fe/H] = —1.6 (Laird 
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et al. 119881 . The halo population in the solar neighborhood is 
not purely a relic of a monolithic, "rapid" collapse (Carney 
et al. I1996> . There have been several suggestions of a two- 
component halo, with a flattened component in the inner halo 
and a more spherical outer halo (Sommer-Larsen & Zhen 1990 
Carney et al. lT996l . 

The Galactic thick disk is the kinematically hottest portion 
of the disk of the Galaxy, with a scale height of 1.0 to 1.5 
kpc and rotates with a velocity of about 170 km s _1 (Gilmore, 
Wyse, & Kuijken 1989 1. The thick disk is usually considered 
to be dominated by stars in the range [Fe/H] > -1 (Freeman 
H9571 . peaking at about [Fe/H] = -0.5 (Carney et al. 1989 1. 
Many workers have claimed the existence of a metal weak tail 
of the thick disk component in the range — 1 .6 < [Fe/H] < — 1 
(MFF; Beers & Sommer-Larsen 1995] Chiba & Beers I^UtM 
MFF found a fraction of 72% of the stars in this metallicity 
range in a "metal-weak thick disk" (MWTD), rotating rapidly 
at y iot « 170 km s i. Another large fraction of MWTD was 
found also by Beers & Sommer-Larsen ( 1995}. Their MWTD, 
rotating at V rot « 195 km s~l, accounts for about 60 % of the 
stars in the range - 1 .6 < [Fe/H] < — 1 in the solar neighbor- 
hood, and it possesses an extremely metal-weak tail down to 
[Fe/H] < -2. Chiba & Beers d2000> estimated the fraction of 
MWTD at about 30 % of the metal-poor stars in the abundance 
range —1.7 < [Fe/H] < -1, which is smaller than the frac- 
tion derived by MFF and Beers & Sommer- Larsen dl995l . but 
larger by ~ 10% than the result of Chiba & Yoshii ( 1998 1 using 
solar neighborhood red giants and RR Lyrae stars. 

The investigation of thick-disk and halo kinematics may 
only be applicable to a specific place in the Galaxy and may 
have fine structure of the velocity distribution smoothed out 
by the velocity resolution of the study (Martin & Morrison 
1998). Here, we study the kinematics of solar neighborhood 
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Fig. 1. Color-Magnitude diagram for all identified CLLA stars. 
Hipparcos parallaxes were used to determine My and its stan- 
dard error. The full lines indicate the zero age main sequence 
(ZAMS) and main sequences of the old open clusters M 67 
and NGC 188. The dashed line is the ZAMS shifted upward by 
AMy = 0.8 mag., used to remove the contamination by sub- 
giants and giants. 
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Fig. 2. Hipparcos trigonometric parallaxes versus the photo- 
metric parallaxes of CLLA, for 539 stars (top-left) and for dif- 
ferent metallicity cuts : [Fe/H] < -1.6 (top-right), [Fe/H] > -1 
(center-left), -1.6 < [Fe/H] < — 1 (center-right), and for stars at 
large distances, 7THi p < 25 mas (bottom-left). The full line is a 
linear fit to the data. 



subdwarf stars based on the sample of high proper motion stars 
by CLLA. CLLA have measured photometric parallaxes, radial 
velocities, and metallicities of mainly A to early G stars, many 
late G and some early K stars in the Lowell Proper Motion 
Catalogue. In total the CLLA sample contains 1464 stars. In 
their paper there are listed 1269 stars with kinematical param- 
eters and 1261 stars with metallicity parameters, and there are 
1447 stars with radial velocities in their catalog. The radial ve- 
locity precision of their sample lies in the range of 0.4 to 1.3 
km s _1 . About 15% of their sample are binaries or multiple sys- 
tems. The typical accuracy of the metallicities was estimated to 
be +0.13 dex. 

The photometric parallax of CLLA was replaced in our 
study by using the high precision parallax of Hipparcos cata- 
logue. We used the Astrometric Catalog TYC2+HIP (Wielen et 
al. |2001 1 for the proper motions of stars with Hipparcos par- 
allaxes. This catalogue is derived from a combination of the 
Hipparcos Catalogue with proper motions given in the Tycho-2 
catalogue with direct solutions (Wielen et al. 2001 1 and previ- 
ous earthbound measurements. We still use the high precision 
radial velocities and metallicities of the CLLA catalogue for 
our study. Previous work using Hipparcos subdwarfs was done 
by Reid (19983 and Fuchs, JahreiG and Wielen JT9981 hereafter 
FJW). In this previous work we discussed the kinematical be- 
haviour of the 560 subdwarfs for which improved parallaxes 
and proper motions were obtained by Hipparcos, in relation to 
their metallicities. In the present paper we increase the size of 
the sample considerably by applying a correction to the photo- 



metric CLLA distances determined using stars with Hipparcos 
parallaxes. 

In Sect. 2 we describe the selection of the sample and the 
photometric correction. In the next section we determine the 
kinematical properties, i.e. the 3D velocities with respect to the 
Sun, and discuss the kinematical behavior of the subdwarfs, in 
particular in relation to their metallicities. Our conclusions are 
summarized in Sect. 4. 



2. Data 

The CLLA data set of 1447 stars has been cross-identified with 
the Astrometric Catalogue TYC2+HIP (Wielen et al. 2001 1 and 
we found 545 stars in common. About 700 CLLA stars did 
not appear in the TYC2+HIP catalogue, but were then cross- 
identified with the Tycho-2 Catalogue (H0g et al. 2000 1, and we 
found 259 subdwarfs with Tycho-2 proper motions. The proper 
motion accuracy of Tycho 2, derived from a comparison with 
the Astrographic Catalogue and 143 other ground-based astro- 
metric catalogues, is about 2.5 mas yr _1 . CLLA used Luyten's 
NLTT proper motions for the calculation of the space velocity 
components. These proper motions have typical errors of 20 to 
25 mas yr . Therefore, TYC2+HIP and Tycho-2 proper mo- 
tions provide an enormous improvement in the accuracy of the 
tangential velocities. 
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2.1. Color Magnitude Diagram 

Fig. 1 shows the color magnitude diagram for all 545 identified 
CLLA stars. The absolute magnitudes and their standard er- 
rors are based on the Hipparcos parallaxes and their errors. The 
B — V colors were taken from the Hipparcos catalogue. Some 
stars for which no distance was given were already recognized 
by CLLA as subgiants. However, we can see clearly from the 
CM-diagram in Fig. 1 that there is still contamination by previ- 
ously undetected subgiants and giants. To avoid these, we have 
removed all stars lying above a line in the CM-diagram defined 
by the zero age main sequence of stars with solar metallicity 
shifted upward by AM V =0.8 mag. About 8% contaminating 
stars were eliminated in this way. 

2.2. Test of Photometric Distances 

We determined the overall correction of the photometric dis- 
tance scale of CLLA by analyzing the parallax difference. 
There are 539 CLLA stars which have both photometric and 
trigonometric parallaxes in our sample (CLLA-TYC2+HIP). 
We compared the Hipparcos parallax with the photometric par- 
allax of CLLA (as shown in Fig. 2). The error bars represent 
both Hipparcos and CLLA parallax errors. A typical error in 
the absolute magnitude of CLLA stars of <Tm v = 0.3 mag is 
assumed. We used a least x 1 method applicable when the data 
have errors in both coordinates. The ^-function is chosen ac- 
cording to Press et al.( 1992 1 
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The slope of the regression b derived from our 539 subd- 
warfs, is b — 1.116 + 0.008. For comparison, a sample with 
7THip < 25 mas leads to a larger correction of b — 1 .169 ± 0.028. 
We also tried cuts in the metallicties, with the result that for 
more metal-poor stars a larger correction was needed. For stars 
with [Fe/H] > -1, which are dominated by thick disk stars, 
-1.6 < [Fe/H] < — 1, and extreme metal-poor stars with [Fe/H] 
< -1.6 we find slopes of the regression line b — 1.093 + 0.010, 
b = 1 .324 + 0.042 and b = 1 .394 + 0.043, respectively. JahreiB, 
Fuchs and Wielen ( 119971 and FJW have found similar correc- 
tions using a smaller sample of subdwarfs. The data points at 
the upper right corner of the first plot (All Data) and the fourth 
plot (-1.6 < [Fe/H] < -1) represent the star HIP 57939, which 
is the nearest star in our sample with 7THi p = 109 mas. No sig- 
nificant changes in the slopes of both plots (less than lcr) are 
found, if we omit HIP 57939 when calculating the slopes. 

All these corrections are used to calibrate CLLA photo- 
metric parallaxes in our sample with no Hipparcos parallaxes 
and about 35 stars with low accuracy Hipparcos parallaxes i.e. 
7THip < 5 mas and ^Hip/ov < 3 mas. 

3. Kinematical Properties 

Using the parallaxes supplied by the Hipparcos catalogue, 
proper motions by the TYC2+HIP and Tycho-2 catalogue, and 
radial velocities given in CLLA, the space velocity components 
U, V, and W, which are directed to the Galactic center, direction 
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Fig. 3. Spatial distribution of the samples CLLA-TYC2+HIP 
(a to c) and CLLA-Tycho-2 (d to f) , respectively. X points to- 
wards the Galactic center, Y in the direction of galactic rotation 
and Z towards the Galactic north pole. 
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Fig. 4. Space velocity components (U, V rot , W) versus metal- 
licity [Fe/H] of the samples CLLA-TYC2+HIP (a to c) and 
CLLA-Tycho-2 (d to f) , respectively. 



of galactic rotation, and north galactic pole, respectively, have 
been calculated with respect to the Sun and then reduced to the 
LSR (local standard of rest). For the latter Delhaye's (1965 1 
values +9, +12, +7 km s~' were adopted for Uq, Vq, Wq, re- 
spectively. Finally, the velocity components were transformed 
onto a frame rotating with circular velocity y c i lc = -220 km 
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Fig. 5. Toomre diagram: (U 2 + W 2 ) 1 ^ 2 versus V rot distributions 
of CLLA-TYC2+HIP and CLLA-Tycho2 samples. The solid 
lines represent total velocities of 50, 100, 150 and 200 km s _1 , 
respectively. 



s 1 relative to the LSR, i.e. the expected rest frame of our 



Galaxy (e.g. Wielen 1986). The rotational velocity is defined 
as V rot = V- Van;. 

Figure0]shows the U, V rot , and W velocities of the samples 
CLLA-TYC2+HIP and CLLA-Tycho2 as scatter plots. The U- 
distribution indicates that the present sample was kinematically 
selected. The CLLA Catalog is based on a proper motion cat- 
alog so that stars with small tangential velocities are missing. 
We can see clearly that for small U-values the diagrams are 
sparsely populated. This is also seen in the V-velocities. The 
stars with metallicities [Fe/H] > - 1 lag on the average by about 
40 km s _1 , i.e. are thick disk stars. The old thin disk stars are 
missing (cf. also Fig|5jl. In the W-velocities no kinematical bias 
is visible; it is apparently lost in projection. However, since we 
are mainly interested in the kinematics of the halo stars, this 
bias is of no consequence in the present context. It is evident 
from this figure that metal-poor stars with [Fe/H] < -1 have 
larger random motions compared with metal-rich ones [Fe/H] 
> - 1 . This shows that the kinematic properties change rather 
abruptly at [Fe/H] « -1 to -2, which is probable the transition 
region from halo to disk component (Ryan & Norris, [T991l and 
Chiba & Yoshii, IT998l . 

The mean motion with respect to the LSR and velocity dis- 
persions were calculated for different groups in [Fe/H] . The re- 
sults are presented in Tabled The most metal-deficient stars in 
the samples, more metal poor than [Fe/H] = -1.6, are dom- 
inated by members of the halo population. These stars ex- 
hibit a radially elongated velocity ellipsoid (cr\j, cry, cr w ) = 
(189 ± 13,97 + 7,98 + 7) and (157 ± 12,87 + 7,77 + 6) km 
s _1 and show no net rotation, (Vrot) = 1 ± 13 and 15 + 14 km 
s _1 for the CLLA-TYC2+HIP and CLLA-Tycho2 samples, re- 



spectively, which are in good agreement with RR Lyrae kine- 
matics of Martin & Morrison ( 1998 ) and Layden et al. {1996 1. 
Chiba & Beers ( 2000) found a lower velocity dispersion in the 
U-direction, (cr v , cry, cr w ) = (141 ±11, 106 ± 9, 94 ± 8) km s _I 
from their 1203 non-kinematically selected stars. 

The velocity dispersion components of the sample in the 
more metal-rich abundance ranges decrease as the contribution 
of the thick disk component progressively increases. In partic- 
ular, for [Fe/H] > —1.0 the contribution of the halo component 
is expected to be negligible. Our CLLA+TYC2+HIP sample in 
this metallicity range has velocity dispersions (<t\j, cry, <x w ) = 
(74 ± 2,50 + 1,37 ± 1) with V rot = 176 km s _1 , which is 
in agreement with thick disk samples of Martin & Morrison 
JT998T ) of RR Lyrae stars and Chiba & Beers ( 120001 of solar- 
neighborhood stars. 

Our CLLA-Tycho2 sample is more sparsely distributed 
than the CLLA-TYC2+HIP sample in the metallicity range 
[Fe/H] > -1. To understand this, we note that the CLLA- 
Tycho2 sample was drawn from the CLLA stars that are not in 
the Hipparcos catalogue. This might imply that mainly CLLA 
stars with magnitudes brighter than 10.5 mag fall in our CLLA- 
TYC2+HIP sample and stars with magnitudes fainter than 10.5 
mag are in the CLLA-Tycho2 sample. Stars with fainter appar- 
ent magnitudes are at larger distances and velocities compared 
to the Hipparcos stars. The minimum distances for each sample 
are 17 and 40 pc for CLLA-TYC2+HIP and CLLA-Tycho2, re- 
spectively. Figure where we plot the spatial distributions in 
X, Y and Z shows this clearly. We can find the minimum tan- 
gential velocity using the minimum distances and mean proper 
motions for both samples, using 



T min 



= 4.74 



iooo' 



(2) 



where (fi), d n 
-l 



n , and Vxmin denote mean proper motions in 
mas yr l , minimum distances in parsecs and minimum tangen- 
tial velocities in km s _1 for each subdwarf sample. We found 
for the CLLA-Tycho2 sample Vrmm > 50 km s _1 , which might 
explain why there are comparatively few thick disk stars in this 
sample (cf. Fig|5}- 



3.1. V rot distributions of subdwarf s 

The corresponding V rot -velocity distributions for the samples 
CLLA-TYC2+HIP and CLLA-Tycho2 are shown in Figs|6]and 
respectively. The first group, [Fe/H] > -1 dex, represents 
what are obviously the thick disk stars. The third group, [Fe/H] 
< - 1 .6 dex, consists of extreme metal-poor stars, dominated by 
members of the halo population. The histograms of the sam- 
ples CLLA-TYC2+HIP and CLLA-Tycho2 can be fitted by 
Gaussian distributions. 

The second group, —1.6 < [Fe/H] < — 1 dex, shows a pe- 
culiar kinematics. A Kolmogorov-Smirnov test, which avoids 
binning of the data, shows that the velocity distribution of the 
very metal-poor stars, [Fe/H] < -1.6, in the combined sample 
is statistically different from the velocity distribution of the in- 
termediate population, - 1 .6 < [Fe/H] < -1.0. The maximum 
deviation of the normalized cumulative distribution between 
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Fig. 6. Rotational velocity (V rot ) distributions of sample 
CLLA-TYC2+HIP grouped according to their metallicities. 
The velocities are reduced to the local standard of rest. 



the two groups is D — 0.241 and thus significantly larger than 
the critical value A).05 = 0.196 (Sachs. [T988> . which leads to a 
rejection of the hypothesis of the statistical similarity of the ve- 
locity distributions of the two groups. Similarly we have shown 
that the velocity distribution of the metal-poor stars is symmet- 
ric with respect to V rot = km s" 1 (D = 0.149, D .05 = 0.282), 
whereas the velocity distribution of the intermediate group is 
asymmetric (D = 0.375, Do.os = 0.300). Thus the intermediate 
group seems to represent a different population of halo stars. 

On the other hand, the asymmetric drift ratio (V)/cr}, = 
-0.007 (CLLA-T YC2 +HIP) is very similar to that of the thick 
disk stars, (V)/cr^ = -0.008. We conclude tentatively that 
the stars in the -1.6 < [Fe/H] < —1.0 metallicity range repre- 
sent a population of the dynamically hot metal-weak thick disk 
(MWTD). 
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Fig. 7. The same as Fig. [6] but for sample CLLA-Tycho2. 



4. Summary and Discussion 

We have analyzed the kinematics of 742 nearby metal-poor 
subdwarf stars from the CLLA catalogue. The subdwarf s were 
cross-identified with the TYC2+HIP and Tycho 2 Catalogues 
to find accurate trigonometric parallaxes and proper motions. 
The accurate Hipparcos parallaxes lead to an upward correc- 
tion factor of 1 1% of the photometric distance scale of CLLA, 
and it was used to correct the photometric distances of CLLA- 
Tycho2 stars. 

The present analysis indicates that the solar neighborhood 
subdwarf stars with [Fe/H] <-1.6 show halo kinematics charac- 
terized by a radially elongated velocity ellipsoid and no signif- 
icant rotation. At a metallicity range of [Fe/H] > -1, our sam- 
ples show disklike kinematics. In the metallicity range -1.6 < 
[Fe/H] < - 1 .0 we found a significant number of stars with kine- 
matics not of halo stars but that of a dynamically not-metal- 
weak tail of the thick disk. 
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Chiba & Beers (20001 obtained a fraction of 30 % of low- 
metallicity stars in their nonkinematically selected solar neigh- 
borhood sample with -1.7 < [Fe/H] < -1.0, which is consistent 
with our result of 18%. Chiba & Yoshii ( 119981 analyzed the 
kinematics of red giants and RR Lyrae stars in the solar neigh- 
borhood based on Hipparcos data. They found in both red gi- 
ant and RR Lyrae samples in the range -1.6 < [Fe/H] < -1.0, a 
fraction of ~ 10 % of stars in a population with a mean velocity 
(V (S >disk = 195kms- 1 . 

We must try to understand the implications of a significant 
population of MWTD stars for theories of the formation and 
evolution of the Galaxy. It should be kept in mind that, although 
the MWTD population may contribute a large fraction of the 
local metal-poor stars, the (inner) halo population is probably 
still dominated by the stars with [Fe/H] < -1.6 within a few 
kiloparsecs of the Sun. Furthermore, although we have empha- 
sized the possible importance of the MWTD population, it cer- 
tainly still appears to be a minor constituent of the entire thick 
disk population (Beers et al. 2002 1. 

If there is indeed a significant fraction of thick disk stars 
with metal abundance -1.6 < [Fe/H] < -1, as we have ar- 
gued, this finding may have significance for formation scenar- 
ios of the Galaxy. An interesting scenario for the origin of an 
MWTD component may be the merging of satellite galaxies 
(Searle & Zinn. ll978l . which are then accreted by a thin, fast 
rotating, possibly metal-poor, Galactic disk (Quinn et al. 119931 
Wyse 2001 1. The dynamical heating of the stellar component of 
this disk in connection with the accretion process produces the 
thick disk. The kinematics of the halo depends on the dynamics 
of the merging satellites, whereas the kinematics of the thick 
disk are determined by the heating of the rotating thin disk. 
Based on this merging picture of galaxy formation, one might 
argue that the "shredded satellite" stars retain a kinematic sig- 
nature distinct from the thick disk part that results from the 
heated thin disk. The kinematic trace of the destroyed satellite, 
which is probably the origin of the MWTD stars, would be vis- 
ible in the mean orbital rotational velocity of stars. Based on a 
spectroscopic survey of ~ 2000 F/G stars 0.5 - 5 kpc above the 
Galactic plane Gilmore et al. (2002i determined a mean rota- 
tional velocity lag of the shredded galaxies of ~ 100 km s _1 . 
The actual lag expected from the shredded satellite depends 
predominantly on the initial orbit and the amount of angular 
momentum transport in the merger process and is not initially 
predictable in a specific case (Gilmore etal.|2002 1. 
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